Kep^ph/01051471 
ANL-HEP-PR-0 1-023 
UAHEP013 

Upsilon Decay to a Pair of Bottom Squarks 

Edmond L. Berger a and L. Clavelli fc 

a High Energy Physics Division, Argonne National Laboratory 
Argonne, IL 60439 
email: berger@anl.gov 
b Department of Physics and Astronomy, University of Alabama 
Tuscaloosa, AL 35487 
email: lclavell@bama.ua.edu 
(May 15, 2001) 

Abstract 

We calculate the rate for T decay into a pair of bottom squarks as a function 
of the masses of the bottom squark and the gluino. Data from decays of the 
T states could provide significant new bounds on the existence and masses of 
these super symmetric particles. 
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Introduction. The possibility of a light bottom squark, b, with mass of order 10 GeV or 
less, is examined in several theoretical and experimental 0-[J papers. If the mass 

is less than half the mass of one of the J PC = 1 T states, then the decay T — > bb might 
proceed with sufficient rate for experimental observation or exclusion of a light b. In this 
paper, we compute the expected rate for T decay into a pair of bottom squarks as a function 
of the masses of the bottom squark and the gluino. The mass of the gluino g enters because 
the gluino is exchanged in the decay subprocesses we compute. The bottom squark and the 
gluino are the supersymmetric partners of the bottom quark and gluon, respectively. 

Measurements at e + e~ colliders do not yet constrain significantly the existence of fe's in 
the mass region < 10 GeV If a light bottom squark (pi) is an appropriate mixture 
of left-handed (L) and right-handed (R) bottom squarks, its tree-level coupling to the neu- 
tral gauge boson Z can be made small by adjustment of the L — R mixing, leading to good 
agreement with observables at the Z [0 . Owing to their p-wave coupling to the intermediate 
photon and fractional charge —1/3, bottom squarks make a small contribution to the inclu- 
sive cross section for e + e~ — > hadrons @,^], in comparison to the contributions from quark 
production, and bb resonances are likely to be impossible to extract from backgrounds j3|. 
The angular distribution of hadronic jets produced in e + e~ annihilation can be examined 
in order to bound the contribution of scalar-quark production. Spin- 1/2 quarks and spin-0 
squarks emerge with different distributions, (l±cos 2 #), respectively. The measured angular 
distribution M is consistent with the production of a single pair of charge- 1/3 squarks along 
with five flavors of quark-antiquark pairs 0. The exclusion of a b with mass 3.5 to 4.5 
GeV does not apply since that analysis focuses on decays into leptons, b — > civ and b — > cl. 
A long-lived b or one that decays via baryon-number-violating .R-parity violating couplings 
would evade the CLEO limitation M. 

The possibility that the gluino may be much less massive than most other supersymmetric 
particles is intriguing from different points of view [|lC| -[T3l. The current status of the very 
low mass gluino scenario, further possibilities, and more complete references are found in 



Ref . |TI| . The production of a pair of gluinos with mass in the range 12 to 16 GeV, decaying 
with 100% branching fraction into a bottom quark and a bottom squark in the mass range 
2 to 5.5 GeV, offers a supersymmetry (SUSY) explanation for the large bottom quark cross 
section at hadron colliders, and it is consistent with hadron collider data on the time- 
averaged B°B° mixing parameter x ■ A renormalization group argument 0] suggests that 
the existence of a light b goes hand-in-hand with a comparatively light g. 

To compute T decay into a pair of bottom squarks, we follow the established color 
singlet approach [|14] -[T6| in lowest-order perturbative quantum chromodynamics (QCD). We 
assume that the T is a bb bound state of bottom quarks and that each b carries half of 
the momentum of the T with 2mb = mx- Refinements may be considered, but this simple 
approach is sufficient for the purposes of this first paper on the topic. 

Notation and Normalization. Define the four- momentum of the T to be pr- The po- 
larization vector of the T is denoted e. The J p = 1~ bound state of the b and b quarks is 
represented as 

_L tt^ ACij , (!) 

where the color matrix Co- = expresses the color singlet system of the b with color 
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index % and b with color index k. The coupling strength of the T to bb is specified in terms 
of an overall constant A which is related to the value of the orbital wave function at the 

wn 



origin m momentum space 



width r^r or the hadronic width rj. In leading order, 



This constant may be expressed in terms of the electronic 



T _ 8irA 2 e 2 a 2 m 

Tfly 

where ej = —1/3. The lowest-order perturbative QCD expression for a three-gluon decay 
width is JT7| 

10(tt 2 - 9) a 3 s (fi) 
T T - 9tt a 2 ' 1 ' 

where the strong coupling a s (/i) is evaluated at a momentum scale /i ~ m T . However, 



next-to-leading order contributions are known to be large |18|]. In our work, we prefer to use 
the electronic width as our source of absolute normalization. 

The T decays into a pair of bottom squarks b of mass m\ carrying four-momenta k\ and 
k 2 respectively. This decay proceeds via the two-parton to two-sparton subprocess 

b + b^b + b, (4) 

with exchange of a t-channel gluino g of mass m§. The subprocess is sketched in Fig. 1. The 
subprocess in which the bb pair annihilates through an intermediate gluon into a bb pair is 
absent because the initial state is a color singlet. 

Mixing. The mass eigenstates of the bottom squarks, b\ and b 2 are mixtures of left- 
handed (L) and right-handed (R) bottom squarks, &l and &r. The mixing is expressed 
as 

61 = cos 9~ b b R + sin 9~ b b L , (5) 
b 2 = - sin 8~ b b R + cos 9- b b L . (6) 

In our notation, the lighter mass eigenstate is denoted b\. For the case under consideration, 
the mixing of the bottom squark is determined by the condition that the lighter b coupling 
to the Z boson be small [0], namely sin^ ~ 0.38. 

Matrix Element and Decay Width. The coupling at the three-point vertex in which a b 
quark enters and a b\ squark emerges (the upper vertex in Fig. 1) is 



.,V2T-[cos9 b P R -sme- b P L ], (7) 



where i and j are the color indices of the incident b and final b, respectively, and a labels the 
color of the exchanged gluino. g s is the strong coupling, a s = g 2 /4:7i. Here, Pl = (1 — 7s)/2 
and P R = (1 + 75) /2. At the lower vertex where an antiquark enters and bi emerges, the 
coupling is 



; V2T^[cos9 b P L -sm9 b P R ], (8) 
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where k and I are the color indices of the incident b and final b, respectively. 

Define t = (ki—pr/2) 2 = — (m^ — 4m 2 )/4. The invariant amplitude for T — > b(ki)+b(k 2 ) 

is 



M = ^™ s AF C Tr ( ^ + mr)[cos0iPR _ sindl p L ^ r _ 2 fa + 2m 5 ][cos0gP L - sin^J) 

(9) 

The color factor F c = ^Tg^jfc, where j, k label the colors of the final squarks. The relatively 
large mass of the exchanged gluino should justify the use of simple perturbation theory to 
compute the decay amplitude. Evaluating the trace and using e • pr = 0, we obtain 

t-m§ 

There is no explicit dependence on the mixing angle. 

After squaring M, averaging over the spins of the T, and summing over the colors of the 
final bottom squarks, we derive 

2 



^,,,,9 32 • 167r 2 , tm T 



At threshold, t -> and |M| -> 0. 

To obtain the decay width of the T into bb, we use the expression 



■>7r ' ' m\ ' 



where \k\ = \\fm\ — Am\. We observe that dT oc |/c| 2 ^ +1 , with i = 1, as expected for the 
p-wave decay of the J p = 1~ T into a pair of scalars. 

Predictions. Our results are presented in terms of a ratio 

nT m T (m 2 Y -4mg) 3 / 2 
^ T« 3 l a em J (t-m 2 ) 2 ' Udj 

so that the absolute value of the normalization constant A cancels. We use the value a s (fi) = 
0.22 ± 0.02 determined from data on T decays and evaluated at the scale /i = m;, = 4.75 
GeV [T9|] . This value of a s (n) is consistent with determinations from other processes at other 
physical scales when all are plotted in common with the predicted 0(a\) QCD dependence 
on ft, with a s (M z ) = 0.1184 and = 213 MeV 

The choice of a value of a s smaller than 0.22 might be deemed more appropriate. 
The presence of a light gluino slows the running of a s (fi). Extrapolation from a s (Mz) 
to a s {rrib) with inclusion of a gluino of mass 16 GeV reduces a s (rrib) to ~ 0.189. The value 



a s (9.45GeV) = 0.163±0.002±0.014 is extracted from data on radiative decays of the T [20 



We comment below on the quantitative changes in our results from use of a s = 0.16 rather 
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than a s = 0.22. A smaller value of a s leads to a more generous range of possible values of 
m g and m~ b . 

As long as mj > \t\ = (m\ — 4ml)/ 4, as i s true f° r the parameter space of interest to 
us, the functional dependence of Eq. (|i~3"D shows that the ratio Rj~ increases with m^. For 
a given bottom squark mass, larger values of Bl= correspond to smaller values of the gluino 

bb 

mass. 

The values of T e i(nS) listed in the Particle Data Group's |HJ] compilation are 1.32 ± 
0.04 ± 0.03 keV, 0.52 ± 0.03 ± 0.01 keV, ~ 0.47 keV, and 0.248 ± 0.031 keV, for n = 1 - 4, 
respectively. The corresponding full widths T T are 52.5 ± 1.8 keV, 44 ± 7 keV, 26.3 ± 3.5 

keV, and 14 ± 5 MeV. A ratio i?X = 1 would fall within a one standard deviation (a) change 

' bb v y b 

of T T in all cases, and a ratio RJ= = 10 is close to or within la except for the T(IS') where 
it corresponds to 25% of r T (lS'). As a working hypothesis, we take R?= = 1 to 10 as a 
definition of the range in which to find solutions for the bottom squark and gluino masses. 
In Fig. 2, we present our results for Rj= as a function of mi and m^. We show sets of 

bb on 

curves for the choices RJ~ = 1 and 10. Within each set, the curves correspond to T(nS) 
with n — 1 — 3, with n increasing from the lowest to uppermost curve. The lines can be 
read as predictions of the width for the corresponding values of bottom squark and gluino 
masses, or as lower limits on the sparticle masses given a known bound on i??=. 

bb 

The requirement RT= < 10 can be satisfied readily within the range of values of the gluino 

bb 

and bottom squark masses favored in Ref. [0. For example, with i?X < 10, m g < 16 GeV is 
achieved for m~ h > 3.5, 3.9, and4.2 GeV at the T(nS) with n = 1 — 3. The more restrictive 
condition RT~ < 1 at the T(1S), along with m g < 16 GeV, would demand that m~ b > 4.5 
GeV, still within the range of acceptable solutions in Ref. [0. While these conclusions are 
consistent with the range of proposed masses in Ref. our analysis of T decays shows 
nevertheless that bounds on the ratio i?X are potentially powerful for the establishment 
of lower bounds on masses. For example, as is evident in Fig. 2, at the T(1S) the choice 
m~ b = 3.0 GeV would require m g > 18 GeV or > 32 GeV if R T = < 10 or < 1, respectively. 

Rather than focusing on the value of i?X we may examine instead the implications of 

bb 

the additional SUSY contribution for the full width of the T. Because full width of the 
T(IS') is the one known most precisely, the most sensitive results on the possible exclusion 
of decays into a pair of 6's with low mass come from measurements at the T(IS'), as long as 
m x(i5) > 2 m b- The hadronic width of the T is calculated in conventional QCD perturbation 
theory [|17||l"8f from the three-gluon decay subprocess, T — > 3g, and, as indicated in Eq. @, 
T 3g oc aj. The SUSY subprocess of interest here adds a new term to the hadronic width 

from T — > b + b. If this new subprocess is present but ignored in the analysis of the hadronic 
width, the true value of a s (n) will be smaller than that extracted from a standard QCD fit by 
the factor (1 — TsusY/rsg)^- For a contribution from the bb final state that is 25% of r T (15') 
{i.e. Rj^ = 10), the value of a s extracted via Eq. (|3]) will be reduced by a factor of 0.9, at 
the lower edge of the approximately 10% uncertainty band on the value [|19| a s = 0.22 ±0.02 
used in our calculation. As remarked above, an even smaller value is extracted from data 
on radiative decays of the T [2Cfl. From the perspective of this examination of T T , we 
deduce that a contribution from bb decays as large as 25% of r T (15') cannot be excluded. A 
thorough analysis would require the computation of next-to-leading order contributions in 
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SUSY-QCD to both the 3g and bb amplitudes and the appropriate evolution of a s (/x) with 
inclusion of a light gluino and a light bottom squark. 

The data sample is largest at the T(4S). However, for a fixed m g = 14 GeV, the 
branching fraction into a pair of bottom squarks is about 10~ 3 for m~ b = 2.5 GeV, and it 
drops to about 10~ 4 for m~ b = 4.85 GeV, below current sensitivity. 

Discussion and Conclusions. Our study of T decay to a pair of bottom squarks demon- 
strates that experimental limits on the ratio R?= could establish significant new lower bounds 

bb 

on the masses of the gluino and bottom squark if m b < my/2. In Ref. JIJ, a supersymmetry- 
based contribution is suggested as an explanation for the observed size of the bottom quark 
cross section and time-averaged B°B° mixing probability observed at hadron colliders. The 
proposal requires a bottom squark in the mass range 2.5 to 5 GeV together with a gluino 
in the range 12 to 16 GeV. In this paper, we show that measurements of RJ= could further 
constrain these allowed regions markedly. For example, with m g < 16 GeV, the very restric- 
tive experimental bound RJ~ < 1 at the T(15) requires m b > 4.5 GeV when a s = 0.22 is 
used in the computation (and m b > 4.35 GeV when a s = 0.16 is used). With m g < 16 GeV, 
the bound RJ= < 10 at the T(nS) would require m b > 3.5,3.9, and 4.2 GeV, respectively, 
for n = 1 — 3 and a s = 0.22 (and m b > 2.6,3.2, and 3.5 GeV for a s = 0.16). The bound 
RJ= < 10 implies that the contribution from the bb final state is as much as 25% of r T (lS'), 

bb 

consistent, within uncertainties, with the determination of a s from T decays. All of these 
values of m b are acceptable from the perspective of Ref. |]J but tend to be towards the upper 
end of the proposed range. Gluinos of mass 12 GeV or less would be incompatible with the 
restriction Rl < 10 unless mg > 4.2,4.55, and 4.75 GeV at the T(nS), n = 1 -3 (3.9, 4.25, 

and 4.45 GeV for a s = 0.16). Gluinos of mass > 35 GeV, or so, are compatible (R?= < 1) 
with current information on the total width of the T states for all m~ b > 2.5 GeV. 

The presence of a SUSY contribution to the total width T T means that the value of 
a s extracted from data will be reduced by the factor (1 — rsusY/r3 5 )3. With a smaller 
a s , the contribution from the conventional 3g decay is reduced, compensated by the SUSY 
contribution, so as to maintain T T . As remarked, use of a smaller value of a s (/x) in our 
calculation allows for a broader range of allowed values (i.e. smaller values) of m g and mg. 
Evolution of a s (fi) from /i = m& to // = mz is slowed by the presence of a light gluino. Thus, 
consistency is achieved because the smaller a s (m&) that allows for a greater contribution from 
T —>■ bb can still be made compatible with the same a s (mz). 

A dedicated experimental investigation of the possible decay T — > bb is warranted. Direct 
observation of the decay requires an understanding of the ways that the 6's may manifest 
themselves 0. If the b is relatively stable, the b could pick up a light u or d and become a 
B~ or B° "mesino" with J = 1/2, the superpartner of the B meson. The mesino has baryon 
number zero but acts like a heavy p - perhaps detectable with a time-of-flight apparatus 
or via dE/dx measurements. R-parity- violating and lepton-number- violating decay of the 
b into at least one lepton, disfavored by CLEO data 0, would imply final states with soft 

leptons. Possible baryon- number- violating R-parity- violating decays of the b are b — ► u + s; 
— > c + d; and — > c + s. These final states of four light quarks might be distinguishable from 
conventional hadronic final states mediated by the three-gluon intermediate state. Studies of 
thrust and other event-shape variables or of exclusive final states may provide discrimination. 
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T's produced in e + e annihilation are polarized transversely because of their coupling to 
the intermediate virtual photon. Correspondingly, the final b squarks from T — ► bb will have 
an angular distribution proportional to sin 2 6. 

A theoretical study similar to that in this paper may be carried out for decays of other 
bottomonium states into pairs of bottom squarks. 
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FIGURES 




FIG. 1. Diagram of the decay of a bound bb system into a bb final state, mediated by single 
gluino exchange. 
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FIG. 2. Loci in the - m g plane for which the rate for decay into a pair of bottom squarks 
is either 10 times T $ for each T(nS), with n = 1 — 3 (lower set of curves) or equal to Y $ (upper 
set). Above the curves, the rate would be less. Within each set of curves, the order is (bottom to 
top) IS,2S,3S. 
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